Magnetic relaxation in a classical spin chain as model for nanowires 
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With decreasing particle size, different mechanisms dominate the thermally activated magneti- 
zation reversal in ferromagnetic particles. We investigate some of these mechanisms for the case 
of elongated, single-domain nanoparticles which we describe by a classical Heisenberg spin chain 
driven by an external magnetic field. For sufficiently small system size the magnetic moments ro- 
tate coherently. With increasing size a crossover to a reversal due to soliton-antisoliton nucleation 
sets in. For even larger systems many of these soliton-antisoliton pairs nucleate at the same time. 
These effects give rise to a complex size dependence of the energy barriers and characteristic time 
scales of the relaxation. We study these quantities using Monte Carlo simulations as well as a di- 
rect integration of the Landau-Lifshitz-Gilbert equation of motion with Langevin dynamics and we 
compare our results with asymptotic solutions for the escape rate following from the Fokker-Planck 
equation. Also, we investigate the crossover from coherent rotation to soliton-antisoliton nucleation 
and multi-droplet nucleation, especially its dependence on the system size, the external field and 
the anisotropy of the system. 
Pacs: 75.10.Hk, 75.40.Mg, 75.40.Gb 



I. INTRODUCTION 



Many novel physical effects occur in connection with 
the decreasing size of the systems which are under in- 
vestigation. Magnetic materials are now controllable in 
the nanometer regime and there is a broad interest in 
the understanding of the magnetism of small magnetic 
structures and particles due to the broad variety of in- 
dustrial applications 0]. Magnetic particles which are 
small enough to be single-domain are proposed to have 
good qualities for magnetic recording and arrays of iso- 
lated, nanometer-sized particles are thought to enhance 
the density of magnetic storage. But on the other hand 
there is an ultimate limit for the density of magnetic 
storage which is given by that size of the particles below 
which superparamagnetism sets in Therefore, the 
role of thermal activation for the stability of the mag- 
netization in nanometer-sized structures and particles is 
studied at present experimentally as well as theoretically. 

Wernsdorfer et al. measured the switching times in 
isolated nanometer-sized particles [^,^, and wires 
For small enough particles they found agreement with 
the theory of Neel and Brown |8| who described the 
magnetization switching in Stoner-Wohlfarth particles by 
thermal activation over a single energy barrier following 
from coherent rotation of the magnetic moments of the 
particle. For larger particles jij and also for wires [^|J|] 
nucleation processes and domain wall motion were found 
to become relevant. 

Asymptotic formulae for the escape rates following 
from corresponding Fokker-Planck equations have been 
derived for ensembles of isolated Stoner-Wohlfarth par- 
ticles as well as for a one-dimensional model 

Most of the numerical studies of the magnetization 
switching base on Monte Carlo methods. Here, mainly 



nucleation phenomena in Ising models have been stud- 
ied 211 but also vector spin models have been used 
to investigate the magnetization reversal in 
systems with many, continuous degrees of freedom qual- 
itatively. However, Monte Carlo methods - even though 
well established in the context of equilibrium thermody- 
namics - do not allow for a quantitative interpretation 
of the results in terms of a realistic dynamics. Only re- 
cently, a Monte Carlo method with a quantified time step 
was introduced Here, the interpretation of a Monte 
Carlo step as a realistic time interval was achieved by a 
comparison of one step of the Monte Carlo process with 
a time interval of a Langevin equation, i.e., a stochastic 
Landau-Lifshitz-Gilbert equation. 

Numerical methods for the direct integration of a 
Langevin equation 1 25 2^ are more time consuming than 
a Monte Carlo method but nevertheless highly desirable 
since here, naturally, a realistic time is introduced by 
the equation of motion. The validity of different inte- 
gration schemes is still under discussion (see for a 
discussion of the validity of Ito and Stratonovich integra- 
tion schemes) and, hence, here, as well as for the Monte 
Carlo methods, the investigation of analytically solvable 
models as test tools for the evaluation of the numerical 
techniques are desirable. 

In this paper we will consider a chain of classical mag- 
netic moments — a system which can be interpreted as 
a simplified model for ferromagnetic nanowires or ex- 
tremely elongated nanoparticles This model is very 
useful for two reasons: i) since it was treated analyti- 
cally asymptotic results for the energy barriers as well as 
for the escape rates are available. Hence, the model can 
be used as a test tool for numerical techniques, ii) De- 
pending on the system size, given anisotropics and the 
strength of the magnetic field different reversal mecha- 
nisms may appear and can be investigated. For small sys- 
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tern sizes the magnetic moments are expected to rotate 
coherently 0| , while for sufficiently large system sizes the 
so-called soliton-antisoliton nucleation is proposed [ p^ . 
Therefore, a systematic nmnerical investigation of the 
relevant energy barriers, time scales, and of the crossover 
from coherent rotation to nucleation is possible. 

The outline of the paper is as follows. In Sec. || first 
we introduce the model and we compare the two differ- 
ent numerical techniques mentioned above, namely the 
Monte Carlo method and the numerical solution of the 
Langevin equation, both of which we will use throughout 
the paper. In Sec. Ill we compare our numerical results 



with the theoretical considerations of Braun, concern- 
ing the energy barriers as well as the mean ffist passage 
times for magnetization switching by coherent rotation 
as well as by soliton-antisoliton nucleation. For higher 
temperatures or driving fields, we also find a crossover to 
multidroplet nucleation, similar to what is known from 
Ising models. In section |^ we summarize our results and 
relate them to experimental work. 



II. MODEL AND METHODS 
A. Model 

Our intention is to compare our numerical results with 
Braun's analytical work which bases on a continuum 
model for a magnetic nanowire. For our numerical in- 
vestigations we use a discretized version of this model 
namely a one dimensional classical Heisenberg model. 
We consider a chain of magnetic moments of length L 
(number of spins) with periodical boundary conditions 
defined by the Hamiltonian, 



x\2 



(1) 



where the Si = fJ-i/^s are three dimensional magnetic 
moments of unit length. The first sum which represents 
the exchange of the magnetic moments is over nearest 
neighbor interactions with the exchange coupling con- 
stant J. The second and third sum represent uniax- 
ial anisotropics of the system where the x-axis is the 
easy axis and the z-axis the hard axis of the system 
(anisotropy constants — O.IJ, = J). These 
anisotropy terms may contain contributions from shape 
anisotropy as well as crystalline anisotropics |]l5|| . Even 
though an exact treatment of dipolar interactions would 
be desirable we let this problem for future work so that 
our results presented here are comparable to the analyt- 
ical work of Braun. The last sum is the coupling of the 
moments to an external magnetic field, where B is the 
induction. 



All our simulations using the two different methods 
above start with a configuration where all magnetic mo- 
ments point into the x direction, antiparallel to the ex- 
ternal magnetic field B = — _Bj;X. 

The time r when the x-component of the magneti- 
zation changes its sign averaged over many simulation 
runs (100-1000, depending on system size and simula- 
tion method) is the most important quantity which we 
determine. In the case where the temperatures is low 
compared to the energy barrier the system is in the 
metastable initial state for a very long time while the time 
needed for the magnetization reversal itself is extremely 
short. In this case t can be considered to be the so-called 
mean first passage time required to overcome the energy 
barrier. For low enough temperatures the mean first pas- 
sage time should also be comparable to the reciprocal 
of the escape rate A which has been calculated for the 
model considered here asymptotically from the Fokker- 
Planck equation in certain limits (for details see |0). 
In the classical nucleation theory the so-called lifetime 
or nucleation time is the time required by the system to 
build a supercritical droplet which from then on will grow 
systematically and reverse the system. The metastable 
lifetime was measured numerically in Ising models by use 
of similar methods ||l6| |l[. For low enough temperature 
T all the different times mentioned above are expected 
to coincide. 



B. Langevin dynamics 



The equation describing the dynamics of a system of 
magnetic moments is the Landau-Lifshitz-Gilbert (LLG) 
equation of motion with Langevin dynamics. This equa- 
tion has the form 



.S,x (S.x (C(i)-||)), 



(2) 



with the gyromagnetic ratio 7 = 1.76 x 10^^(Ts)~^, 
the dimensionless damping constant a. The first part 
of Eq. H describes the spin precession while the second 
part includes the relaxation of the moments. In both 
parts of this equation the thermal noise dit) is included 
representing thermal fluctuations, with (dit)) = and 
(C- (t)C' (i')> = 26^,SkiSit ~ t')2akBT^is/l where i,j de- 
note the lattice sites and fc, I the cartesian components. 
Eq. H is solved numerically using the Heun method which 
corresponds to the Stratonovich discretization scheme 
[p6| . Note, that since the Langevin dynamics simulations 
are much more time consuming than the Monte Carlo 
simulation we use this method here mainly for compar- 
ison and we will present only data for the relaxation on 
shorter time scales. 
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C. Monte Carlo simulations 



III. RESULTS 



It is inconceivable that the Langevin dynamics simu- 
lation can be used over the whole time scale of physical 
interest so that we simulate the system also using Monte 
Carlo methods ||2^ with a heat-bath algorithm. Monte 
Carlo approaches in general have no physical time asso- 
ciated with each step, so that an unquantified dynamic 
behavior is represented. However, we use a new time 
quantified Monte Carlo method which was proposed in 
|P4| where the interpretation of a Monte Carlo step as a 
realistic time interval was achieved by a comparison of 
one step of the Monte Carlo process with a time interval 
of the LLC equation in the high damping limit. We will 
use this algorithm in the following. The trial step of this 
algorithm is a random movement of a magnetic moment 
within a cone of size r with 



We investigate the influence of the system size on the 
occurring reversal mechanisms. There are two extreme 
cases of reversal mechanisms which might occur in our 
model namely coherent (or uniform) rotation and nucle- 
ation. For small system sizes all magnetic moments of 
the particle rotate uniformly in order to minimize the 
exchange energy. For larger system sizes it is favorable 
for the system to divide into parts of opposite directions 
of magnetization parallel to the easy axis which mini- 
mizes the anisotropy energy. This is a magnetization re- 
versal driven by nucleation and subsequent domain wall 
propagation. The crossover between these mechanisms is 
discussed later in this section. First, we study the differ- 
ent reversal mechanism and compare our numerical data 
with theoretical formulae. 



^2 ^ 20kBTa^ 



(1 



(3) 



Using this algorithm one Monte Carlo step represents a 
given time interval At of the LLC equation in the high 
damping limit as long as At is chosen appropriately (for 
details see [pi). 
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FIG. 1. Reduced mean first passage time t^/^s vs. damp- 
ing constant a. The data are from Monte Carlo and Langevin 
dynamics simulations for ksT = 0.025J and L = 80. Solid 
lines are guides to the eye. 

To test the algorithm, in Fig. |l|the a dependence of the 
mean first passage time of the Monte Carlo data is com- 
pared with the data of the Langevin dynamics simulation. 
Each data point is an average over 1000 independent 
runs. For low values of the damping constant the data do 
not coincide while in the high damping limit Monte Carlo 
and Langevin dynamics data converge. Hence, through- 
out this work we will use a = 4, which is large enough 
so that our Monte Carlo simulation and the numerical 
solution of Eq. ^ yield identical time scales. Even though 
this is an unphysically large value for a we can compare 
our results with the analytically obtained high damping 
asymptotes. 



A. Coherent Rotation 

In the case of small chain length the magnetic moments 
rotate coherently while overcoming the energy barrier 
which is due to the anisotropy of the system. The first 
theoretical description of the coherent rotation of elon- 
gated single-domain particles was developed by Stoner 
and Wohlfarth [E8[ . The belonging energy barrier is 



AE„ = Ld4l-hy' 



(4) 



with h = fisBx/{2dx)- Neel expanded this model for the 
case of thermal activation and Brown j|] calculated 
the escape time 



. exp 



kBT 



(5) 



following a thermal activation law. The prefactor r* is 
not a simple constant attempt frequency as claimed fre- 
quently by several authors but a complicated function 
which in general may depend on system size, tempera- 
ture, field and anisotropics. For our model it is 



7r(l+a^) 



a{djl-h^)-d,)+y/(a{d^{l-h^)+d,)y+id^d,{l-h^) 



(6) 



as was calculated from the Fokker-Planck equation . 

The main assumption underlying all results above is 
that the system can be described by a single degree of 
freedom, namely the magnetic moment of the particle 
which must have a constant absolute value. Then, the 
equations above should hold for low enough temperatures 
k^T ^ AE„ and for fisB^ < Sd^,. For larger fields the 
energy barrier is zero, so that the reversal is spontaneous 
without thermal activation (non-thermal magnetization 
reversal) . 

In the following, our numerical results for the mean 
first passage times t for coherent rotation are compared 
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with the equations above. Fig. |^ shows the temperature 
dependence of r for a given value of the external magnetic 
field h = 0.75 and three different system sizes. For tem- 
peratures k^T < Ai^cr our data confirm the asymptotic 
solutions above for smaller system sizes. 
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FIG. 2. Reduced mean first passage time r-y / fj,s vs. inverse 
temperature. The data are from Monte Carlo (open symbols) 
and Langevin dynamics (filled symbols) simulations for differ- 
ent system sizes. The solid lines are the asymptotic formulae 
for the escape times (see Eq. ^ and Eq. ^ . 

For the largest system shown here {L — 16) the nu- 
merical data are systematically lower than the theoreti- 
cal prediction. Obviously, the prefactor r*^ depends on 
the system size in contradiction to Eq. ^ while the en- 
ergy barrier (Eq. ^) is obviously still correct (the slope of 
the data). This size dependence of r*j. can be explained 
with the temperature dependence of the absolute value 
of the magnetic moment of the extended system. In an 
extended system, the energy barrier for vanishing field, 
dxL, is reduced to dxC^iiSf)^) which expanded to first 
order of the temperature can lead to corrections of the 
form d^CEriSf)^) ~ dxL{l-aT). Including this in Eq. | 
remarkably leads to an effective reduction of the prefactor 
T*j. by a factor exp(— ad^^L) and not to a reduction of the 
effective energy barrier AEd-. To conclude, in extended 
systems with many degrees of freedom the reduction of 
the magnetization leads effectively to a size dependent 
reduction of the prefactor of the thermal activation law 
for low temperatures. 



B. Soliton-Antisoliton Nucleation 

With increasing length of the chain a different rever- 
sal mechanism becomes energetically favorable, namely 
the so-called soliton-antisoliton nucleation proposed by 
Braun [Q. Here, during the reversal the system splits 
into two parts with opposite directions of magnetization 
parallel to the easy axis. These two parts are separated 
by two domain walls with opposite directions of rotation 
in the easy x — y-plane (a soliton-antisoliton pair) which 
pass the system during the reversal. 



Note, that we consider periodic boundary conditions, 
otherwise also the nucleation of one single soliton at one 
end of the chain would be possible. The energy barrier 
Aii^nu which has to be overcome during this nucleation 
process is [Ol 



AEnu ^ V^Jdx (4tanhi? - ihR) , 



(7) 



with R = arcosh( -^/T/Ti) . In the limit of small mag- 
netic fields, h 0^ this energy barrier has the form 
Aii^nu — 4:y/2Jdx which represents the well-known en- 
ergy of two domain walls. The corresponding mean first 
passage time follows once again a thermal activation law 
(see Eq. H). Since our Monte Carlo simulations are for 
rather large damping, a = 4, we compare our numeri- 
cal data with the prefactor obtained in the over-damped 
limit (Eq. 5.4 in |l4|) which in our units is 



r.t„7_ ^3/2(1 -fa^)(fcBr)V^(2J)V4 
f^s 16a L dV'^\EoiR)\tanhR^/^smhR 



(8) 



The eigenvalue Eq{R) has been calculated numerically in 
In the limit ft. ^ 1 it is |£;o(-R)| ~ 3i?2. The 1/L de- 
pendence of the prefactor reflects the size dependence of 
the probability of nucleation. The larger the system the 
more probable is the nucleation process. Furthermore, 
the prefactor has a remarkable VT dependence leading 
to a slight curvature in the semi-logarithmic plot of the 
thermal activation law (see Eq. |5|) . 




FIG. 3. Reduced mean first passage time tj/ jj,s vs. inverse 
temperature AEnu/ksT for two different system sizes. The 
data are from Monte Carlo (open symbols) and Langevin dy- 
namics (filled symbols) simulations. The solid lines are the 
asymptotic formulae for the inverse escape rates for a single 
soliton-antisoliton nucleation process (see Eqs. |5|, ^, and g) 
and multiple nucleation (see Eq. pj] ). 

Fig. H shows the temperature dependence of the re- 
duced mean first passage time for h = 0.75 and two dif- 
ferent system sizes. The formulae above are confirmed 
only for rather low temperatures (A£'„u/fcBT > 8 for the 
smaller system and AE^^/kBT > 10 for the larger sys- 
tem). In the range of intermediate temperatures (but 
still k^T < AE^u) the numerical data of both, Langevin 
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dynamics and Monte Carlo simulations deviate from the 
formulae above. Interestingly, in this region the mean 
first passage times do not depend on system size in con- 
trast to Eq. ^ where a 1/ i-dependence occurs due to the 
size dependence of the nuclcation probability. 
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FIG. 4. x-component of the magnetic moments of the chain 
at the time r after single soliton-antisoliton nucleation (up- 
per diagram, L = 80, fceT = 0.0019J, h = 0.95) and multiple 
nucleation (L = 120, ksT = 0.038J, h = 0.95) 

In order to understand this effect occurring in the 
intermediate temperature range. Fig. ^ shows the a;- 
component of the magnetic moments of the chain at the 
mean first passage time r. The upper diagram shows 
the pure soliton-antisoliton nucleation occurring at low 
temperatures. Due to thermal fluctuations a nucleus is 
originated and two domain walls pass the system. The 
lower diagram shows that alternatively also several nuclei 
may grow simultaneously. Obviously, depending on tem- 
perature (and also other quantities, like system size and 
field) with a certain probability many nuclei may arise at 
the same time. This is a multiple nucleation which was 
investigated mainly in the context of Ising models ]29| ] 
where it is called multidroplet nucleation. 



C. Multidroplet Nucleation 



where v is the domain wall velocity and D the dimension 
of the system. Hence, the time when half of the system 
is reversed is given by 

For the one dimensional system which we consider in this 
paper the lifetime is given by 



Lt,* 



■ exp 



A^nu 

2kBT' 



(11) 



This means that the (effective) energy barrier for the mul- 
tidroplet nucleation is reduced by a factor 1/2 and the 
prefactor does no longer depend of the system size 
since t*^ for the soliton-antisoliton nucleation has a 1/L- 
dependence (see Eq. |8|). Since we do not know the value 
of the domain wall velocity we had to fit this parameter 
which influences the prefactor for the comparison in 
Fig. 1^ in the intermediate temperature range. Neverthe- 
less, in Fig. H the reduction of the slope of the curve by 
the factor 1/2 in the multidroplet regime is confirmed. 
Also, it can be seen, that the prefactor is not dependent 
on the system size since the data for both sizes coincide 
as is explained by the considerations above (we used the 
same value for v of course for both curves, since v should 
not depend on the system size). 

A similar crossover from single to multidroplet exci- 
tations was observed in Ising models, field dependent 
|p7|-p^ as well as temperature dependent Compar- 
ing the mean first passage times for single and multiple 
nucleation we get for the intersection of these two times 
the crossover condition 



Lsni = V VT*^Lsm exp 



2fcBT- 



(12) 



The corresponding time is Lgm/v — the time that a do- 
main wall needs to cross the system. This results is also 
comparable to calculations in Ising models [0 . 



The mean first passage time for the multidroplet 
nucleation can be calculated with the aid of the classical 
nucleation theory [^^. In the case of multiple nucle- 
ation many nuclei with a critical size originate within 
the same time. These supercritical nuclei grow and join 
each other leading to the magnetization reversal. The 
mean first passage time for this process is determined by 
the probability for the occurrence of supercritical nuclei 
1 /tuu and the time that is needed for the growth of the 
nuclei. Let us assume that the radius of a supercritical 
nucleus grows linearly with the time t. Then Tmn is given 
by the condition that the change of the magnetization 
AM equals the system size p9|. 



AM(rir 



{2vty 



■dt = L 



D 



(9) 



D. Crossover 

In this section we are interested in the crossover be- 
tween coherent rotation and soliton-antisoliton nucle- 
ation. Therefore, we compare the energy barrier of 
soliton-antisoliton nucleation (see Eq. ^ with the energy 
barrier of coherent rotation (see Eq. |4|) in order to obtain 
a condition for the mechanism with the lowest activation 
energy. The resulting crossover line Lc has the form 



y2Jd;(4tanhi? - AhR) 
rf.(l-/i)2 ■ 



(13) 



For vanishing magnetic field it is Lc = 4y and we get 
the simple condition that four domain walls (not two!) 
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have to fit into the system. Note, that Braun also de- 
termined this crossover line from a slightly different 
condition. He calculated the system size below which no 
nonuniform solution for the Euler-Lagrange equations ex- 
ists and he obtained a limit which is similar to Eq. ^ but 
a little bit lower. 

A diagram showing which reversal mechanisms occur 
in our model depending on system size and field is pre- 
sented in Fig. ||. The crossover line Lc derived above sep- 
arates the coherent rotation region from that of soliton- 
antisoliton nucleation. For h > 1 the reversal is non- 
thermal. In the nucleation region, for larger fields a tem- 
perature dependent crossover to multiple nucleation sets 
in. However, the lower the temperature the more van- 
ishes this region. As an example we show in the figure 
the line Lgra which separates single from multiple nucle- 
ation for k^T — 0.006J and under the assumptions that 
the domain wall velocity is proportional to the field [0 . 
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FIG. 5. The diagram shows the regions of different rever- 
sal mechanisms. The lines corresponds to Eq. |l| and 
The points are data from Monte Carlo simulations confirm- 
ing Eq. |l|. 

In order to confirm Eq. |l^ numerically we determine 
the mean correlations of the y-components of two spins 
located in a distance of L/2 (the maximum distance in a 
system with periodic boundary conditions) at the time r 
during the reversal 

1 ^ 

c,(L/2)=[-^5r(r)5,f^^/,(r)]^^. 

i=l '^^ 

This quantity characterizes the reversal mechanism. 
Fig. ^ shows the system size dependence of Cy{L/2) for 
different values of the magnetic field. For coherent rota- 
tion the limiting value of Cy{Ll2) is 1 since at the time r 
all moments point into the direction (small systems). 
On the other hand, for nucleation the system is mainly 
split into two parts where the moments point into ±a;— 
direction. In this case the correlations in y— direction are 
zero (larger system size). As a criterion for the crossover 
between nucleation and coherent rotation we use a value 
of 1/2 for the correlation function and define the corre- 



sponding system size as the crossover length. This anal- 
ysis leads to the numerical data shown in Fig. ^. 
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FIG. 6. System size dependence of Cy{L/2) for different 
magnetic fields. ksT = 0.1 ... 0.004 J depending on h. 

The MC data confirm the theoretical crossover line ex- 
cept for the largest magnetic field. Here, the numerical 
data differ from the theoretical crossover line which di- 
verges in the limit /i — > 1. In this region a mixed reversal 
mechanism appears: first the magnetic moments rotate 
coherently up to a certain rotation angle and then an in- 
stability sets in and a (restricted) soliton-antisoliton pair 
arises. Obviously, the coherent rotation is unstable in 
this region. A hint for the existence of this mixed mech- 
anism can also be found in Fig. ||: for the field h = 0.95 
the correlation is not zero for large system sizes but re- 
mains finite due to the fact that the whole system first 
rotates by a certain angle before the nucleation sets in. 



IV. CONCLUSIONS 



We investigated the magnetization switching in a clas- 
sical Heisenberg spin chain which we consider as simple 
model for nanowires or elongated nanoparticles and as a 
test tool for numerical techniques since analytical expres- 
sions for the relevant energy barriers and time scales exist 
in several limits. Numerically we used Monte Carlo meth- 
ods as well as Langevin dynamics simulations and we 
confirmed that the Monte Carlo algorithm we use yield 
time quantified data comparable to those of the Langevin 
dynamics simulation in the limit of high damping. 

Varying the system size we observed different reversal 
mechanisms, calculated the mean first passage times for 
the reversal (i.e. the time scale of the relaxation process) 
and compared our results with analytical considerations. 
The comparison of our numerical data with theoretical 
considerations confirms that both of our numerical tech- 
niques can be used to study the relaxation behavior of 
magnetic systems qualitatively and quantitatively. 
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FIG. 7. Reduced mean first passage time T'y/yLs vs. sys- 
tem size L for two different temperatures [k^T = 0.024J and 
k^T — 0.016 J). The solid lines are the appropriate formulae 
and the data are from Monte Carlo simulations. 

Fig. summarizes our results. It shows the system 
size dependence of the reduced mean first passage time 
for two different temperatures. For small system sizes the 
spins rotate coherently. Here, the energy barrier (Eq. |^) 
is proportional to the system size leading to an increase 
of T with system size. Following Eq. |^ the prefactor of 
the thermal activation law should not be L dependent. 
However, here we found slight deviations from the asymp- 
totic expressions stemming from the non-constant mag- 
netization of extended systems. In the regime of soliton- 
antisoliton nucleation the energy barrier (Eq. |^) does not 
depend on the system size but the prefactor (Eq. ||) has a 
1/L dependency. Interestingly, this leads to a decrease of 
the mean first passage time with increasing system size. 
Therefore, there is a maximum relaxation time close to 
where the crossover from coherent rotation to nucleation 
occurs (see Eq. |l^). This decrease ends where the so- 
called multidroplet nucleation sets in (Eq. |lj). From 
here on with increasing system size the mean first pas- 
sage time remains constant (Eq. pi] ). 

Note, that qualitatively the same behavior can be 
found in the particle size dependence of the dynamic co- 
ercivity. Solving the equation describing the thermal ac- 
tivation (Eq. ^ in the three regimes explained above for 
/i(L,T = const), one finds an increase of the coercivity 
in the coherent rotation regime, a decrease in the nucle- 
ation regime, and at the end a constant value for mul- 
tiple nucleation. This is qualitatively in agreement with 
measurements of the size dependence of barium ferrite 
recording particles . 

Comparing the energy barrier for the nucleation pro- 
cess with experimental values (Eq. |^, one should first 
mention that in an experimental system nucleation usu- 
ally will start a the sample ends. This will reduce the 
energy barrier by a factor of 2 since without period- 
ical boundary conditions one does not need a soliton- 
antisoliton pair but just one single excitation. The en- 
ergy barrier for the nucleation process (Eq. [t]) was com- 
pared in with energy barriers measured in isolated 
Ni-nanowires |^,^ finding the experimental energy bar- 



rier reduced by a factor of 3. This agreement is rather 
encouraging taking into account that in realistic systems 
the energy barrier might be reduced depending on the 
form of the sample ends. 

Regarding the crossover from coherent rotation to nu- 
cleation a similar analysis, even though less rigorous, 
has been performed for three dimensional systems [ pO| . 
Here, in the limit h ~* the crossover diameter is 
Lc = 3/2 \J2J jdx instead of 2 ^2J jd^ for the spin 
chain with open boundary conditions. As discussed in 
po| for CoPt particles the crossover length should be of 
the order of 30nm which is also roughly in agreement 
with experiments on Co particles 

We should mention that the model treated here con- 
tains magnetostatic energy only in a local approxima- 
tion ||l5| . Arguments against this treatment have been 
brought forward |3^] so that the question wether an ex- 
plicit inclusion of the dipole-dipole interaction would lead 
to additional reversal modes (e. g. curling) - also depend- 
ing on wether the diameter of the nanowire is larger or 
smaller than the exchange length - remains open. Calcu- 
lations following these lines are therefore highly desirable. 
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